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|5|, |6|] predict similar masses with 
. The Casimir scahng hypothesis 



I. INTRODUCTION 

Lattice calculations [l[ and other models of QCD predict that the Casimir scaling 
hypothesis is well verified in QCD, that is to say that the potential between two opposite 
colour charges in a colour singlet is proportional to the value of the quadratic Casimir 
operator. This mechanism can be tested for instance in two-^luon glueballs. Direct evidence 
for such a system is still controversial, but several mode 
values close to some of the experimental candidates 
can also be tested in another system whose colour-spin structure is similar to the one of a 
glueball. Let us consider an hybrid meson containing a colour octet spin one cc pair and a 
gluon. Due their very heavy masses, the charm quarks can be assumed fixed at the centre 
of mass with the gluon orbiting around. A gluon attached to a heavy colour octet is called 
a gluelump, so the gcc system considered here is called hybrid gluelump as in Ref. j^. 

In our model, a quark and an antiquark are attached by a fundamental string while two 
gluons in a glueball are attached by an adjoint string. We assume here that the three-body 
nature of the confinement implies that the gluon and the colour octet pair are attached by 
two fundamental strings in a hybrid gluelump. Within these hypothesis, we compute the 



masses of glueballs and hybrid gluelumps with a simple effective model of QCD jlOl . Ill |. 
This semirelativistic potential model is described in Sec. |n] It depends on few parameters 
which are fixed on well known mesons in Sec. IIIII Masses of glueballs and hybrid gluelumps 
are computed respectively in Sec. US and El We focus our attention on the first orbital 
excitations, neglecting the details of the spin dependent interactions, assumed weak with 
respect to the confinement. Some concluding remarks are given in Sec. IVII 

II. HAMILTONIAN 

Starting from the QCD theory, a Lagrangian for a system of two confined spinless colour 
sources can be derived taking into account the dynamical degrees of freedom of the string, 



with tension a, joining the two particles, the rotating string mode fRSM) 10|, This 



model is completely equivalent to the relativistic fiux tu 



fields appearing in the RSM are properly eliminated 



equations of these models are difficult to solve jl^ . Ilfil | . So we will use an approximation of 



3e model 12], once the auxiliary 



. The nonlinear coupled 



2 







these effective QCD theories |17|,|18|, which is simply given by the following spinless Salpeter 
Hamiltonian 



Ho = y]? + ml + ^J]? + ml + ar, (1) 
completed by a perturbative part due to the motion of the string 

aL{L + l) 



str 



r 



2/ii/i2 [12/ii/i2 + 4(/ii + fi2)crr + (o-r)^] ' 
The quantity /Xj appearing in the above equation is a kind of constituent particle mass given 
by 



f^i = (^yp^ + mf^ , (3) 

in which the average value is taken for an eigenstate of the Hamiltonian (P). The constituent 
mass is then state-dependent. Equation © is obtained by taking the slow motion limit of 
Eqs. (33) in Ref. Q]. The contribution AMgtr to the mass due to the string can be computed 
with a good precision by the following approximation Q, Q 

2/ii/i2 [12/ii/i2 + 4(/ii + /i2)a(r) + (a(r))2] ' ^ ^ 

It is worth noting that, in this formalism, the particle masses mi and m2 are the current 
ones. So, in the following, we will consider that the mass m„ of the n quark (n stands for u 
or d) and the mass rUg of the gluon are vanishing. 

A significant contribution to hadron masses is also given by the one-gluon exchange 
(OGE) mechanism between colour sources. At the zero order (neglecting the spin), the 
interaction has the following form 

A/7coui = -K— . (5) 
r 



K is a colour factor given by 



f^ = l (Ci2 - Ci - C2) , (6) 



where C12 is the colour Casimir operator for the pair 12 and Ci is the corresponding one for 
the particle i. 
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If Mo is the solution of the eigenequation i^o|0) = ^o|0)) the total mass M for the system 
of the two colour sources is given by 



M = Mo + AM3tr + (Affcoui). (7) 

In this formula, the Coulomb-like contribution is computed as a perturbation. It is shown 
in R.f. Q. that this approximation . very good, especially for states with high angular 

momentum L. All computation are performed using the Lagrange-mesh method |19| . 

Recently, it was shown that the quark self-energy (QSE) contribution, which is created 
by the colour magnetic moment of the quark propagating through the vacuum background 
field, adds a negative constant to the hadron masses 20|. Its negative sign is due to the 
paramagnetic nature of the particular mechanism at work in this case. Other contributions 
due to quark spin (spin-spin, spin-orbit) also exist but they are smaller that the QSE one 



2jJ , and they are neglected in this work. The QSE contribution ASi for a quark of current 



mass nii is given by 

AS. = -/JLMI^iZ^) . (8) 

The rj function is such that r]{0) = 1, and its value decreases monotonically towards with 



increasing quark mass 1^. 6 is the inverse of the gluonic correlation length, and its value 
is estimated about 1.0-1.3 GeV. As the meson masses vary very little with this parameter 
jisl ]. we fix the 6 value at 1 GeV (meson masses do not depend on 6 when m = 0). The 
actor / has been computed by lattice calculations. First quenched calculations gave f = 4: 
2^ . A more recent unquenched work gives / = 3, the value that we choose in this 
work. If a hadron contains q quarks, a contribution J2i=i AS'j must be added to its mass. 
Strong theoretical and phenomenological arguments indicate that gluons do not bring any 



contribution of self energy 



24 



III. PARAMETERS 



As we are mainly interested in describing the main features of spectra, we can use a 
Hamiltonian of spinless particles since spin effects are generally an order of magnitude smaller 
than orbital or radial excitations. Moreover, we will consider hadrons for which the spin- 
dependent part of the Hamiltonian is the weakest possible. 
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TABLE I: Experimental masses in GeV for some light mesons (nn), D* mesons (cn), and charmo- 
niums (cc), with maximal S and J = L + S 30]. The mesons ^(4040) and 7/'(4415) are assigned 
to = 2 and = 3 states respectively. 



N 


= 0, J^(^) l-(-) 2+(+) 


3-{-) 4+{+) 


nn 


7 = 


0.783 1.275 


1.667 2.034 


nn 


/ = 1 


0.775 1.318 


1.689 2.001 


nn 


av. 


0.777 1.307 


1.684 2.009 


cn 




2.010 2.460 




cc 




3.097 3.556 




J = 


= 1, N 


1 


2 3 


cc 




3.097 3.686 4.040 4.415 



The raise of degeneracy between spin and spin 1 mesons can be due to a spin-spin 
potential coming from the OGE or due to an instanton induced interaction 26|, |27|. For 
spin 1 mesons, the spin-dependent interaction is small in the former case or vanishing in the 
latter case. The spin-dependent part of the gluon-gluon interaction has a small contribution 
for glueballs with spin 2 j^. If we consider a hybrid gluelump containing a cc pair with 
spin 1, we minimize the spin effect for the heavy quarks. As such a system is similar to 
a glueball from the point of view of the spin-colour structure, we can also expect that the 
spin contribution will be the smallest if we study spin 2 hybrids. For all these reasons, as in 
Ref. we will only consider systems with maximal value of the total spin S: J = L + 1 
for mesons and J = L + 2 for glueballs and hybrid gluelumps. 

A correct determination of the parameter is crucial to predict new states. For theoretical 
reasons, we have chosen to take m„ = 0, mg = 0, / = 3 and 6 = 1 GeV. Now, we will fix 
the values of the remaining free parameters with well established meson states. As, the light 
nn spin 1 mesons can be isoscalar (/ = 0) or isovector (/ = 1), we will use an average mass 
defined, as usual, by 

M„,(av.) = M..(/ = 0)+3M„„(/=l) 

At this stage, it is interesting to obtain an approximate analytical formula for light nn 
mesons. For such a system, mi = m2 = and /ii = = = (VW)- It can then be shown 



5 



that Mq = 4/i with /i oc ^Ja [14]. Using the semirelativistic virial theorem j28[, we find 
(r) = 1[ijo. With the following approximation (1/r) ^ 1 we obtain finally 

M»«v^f4.-^^^-^5^-i) (10) 

in which = /x/i/a is independent of a. Let us note that an approximate value for v is 
given by [l^ 

..(^)"\ (U) 

where eAr^ is the solution of the dimensionless Hamiltonian + |x| in which q and x are 
conjugate variables. The accuracy of formula (llUj) with the approximation ()lip is about 
5-10% |l7|. 

The contribution of the motion of the string to the meson mass is small, since the ratio 
AMstr/Mo is at most 1/16 118]. Formula ()1()|1 shows that the contribution of the Coulomb 
term and of the self energy shift jthe square mass M'^^ of the meson without modifying the 
slope of the Regge trajectories ll8]. This slopes depend mainly of the string tension which 
we call a for a quark-antiquark bound state. With a fit, it is only possible to determine the 
following combination of parameters — j jn. But since the value of / is fixed, it is 

possible to compute as- 

We fix the parameter by fitting the exact results of Hamiltonian (P) on experimental data 
[3^ (see Table H)). Using the Regge trajectory for light meson, we find a = 0.175 GeV^ and 
as = 0.10. The value obtained for as is small with respect to values found in other potential 
models of meson. For instance, in Ref. ^3]; a value of 0.39 is obtained because the meson 
masses are fitted on the centre of gravity of spin 1 and spin mesons. In this case, the 
average mass of the J = 1 nn state is 0.612 GeV, a value quite different from the one we 
choose, 0.777 GeV. The mass of the quark c is determined by computing the masses of D* 
mesons, which are systems dynamically closer to hybrid gluelumps than charmonium states. 
All parameters are gathered in Table ITTl 

The quality of a fit is estimated by computing for each state the relative error on binding 
energy. The theoretical (experimental) binding energy for a meson is given by the theoretical 
(experimental) mass minus the theoretical current masses of the quark and the antiquark 
(0 for the light nn mesons, rric for the D*, 2mc for the charmoniums) . Quite good meson 
masses are obtained, taking into account the simplicity of the model (see Table ITTT]) . The 
error is around 1% for light nn mesons (see Fig. Q), around 3% for the D* (see Fig. |21); 



TABLE II: Parameters of the model. The values of m^, rUg, nice, f, and 6 are fixed by theoretical 
or phenomenological considerations (see text). 



m„ = 


a = 


0.175 GeV2 


mg = 


as = 


= 0.10 


nic = 


1.300 GeV / = 


3 


rricc = 


: 2 Tile S = 


1 GeV 




FIG. 1: Regge trajectories for light mesons (average between 7 = and 7=1) and glueballs 
(with or without Coulomb interaction). These states are characterized by = 0, S maximal, and 
J = L + S. 

around 10% for the mesons J/ip and Xc2 (see Fig. |21), but less than 5% for radial excitations 
of the J/^ (see Fig.E)). 



IV. GLUEBALLS 



In this paper, we assume that the string tension between two colour sources is controlled 
by the Casimir scaling. So, contrary to what is done in Ref. P|, we assume that the ratio 
between the string tension in a glueball and the string tension in a meson is given by the 
ratio of the Casimir operators. In a two gluon system, we have then cr = 9 a/4, and the 
legge slope is increased by the same factor. Moreover, following the theory of Simonov 
24 . l25l |. there is no self-energy contribution coming from gluon. Let us note that, in the 
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D* (2010) 



A — hybrid (A- = 0) 



12 3 4 

J 

FIG. 2: Masses for cn, cc, and hybrid gluelumps (with or without Coulomb interaction). These 
states are characterized by = 0, S" maximal, and J = L + S. 




FIG. 3: Masses of cc J = 1 mesons as a function of the quantum radial number N. The states 
V'(4040) and ■(/'(4415) are assigned to ip{3S) and Tp{4:S) respectively. 

same framework 10,^^, the constituent mass of a particle is directly proportional to ^/a 
when its current mass m is vanishing. If we compute fig from a two-gluon glueball and fin 
from a light meson, we have fig/^n = 3/2. A gluon cannot decay spontaneously into a nn 
pair. 

At the lowest order, the colour Coulomb factor in such a system is given by the formula ©, 
that is to say k = 3. Nevertheless, in some modelj_it is found that the Coulomb force in 
glueballs could be damped or strongly reduced 2^. So, we have considered both cases 
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TABLE III: Predicted masses in GeV for nn (average between 7 = and 7 = 1), cn, cc, glueballs, 
and hybrid gluelumps. These states are characterized by = 0, S" maximal, and J = L + S. The 
asterisk indicates that the states is calculated with k = 0. 

J 1 2 3 4 5 6 

nn 0.767 1.306 1.670 2.024 2.304 2.555 

Glueball 1.820 2.976 3.793 

Glueball* 1.981 3.059 3.855 

cn 2.033 2.436 2.747 3.010 

cc 3.146 3.449 3.693 3.908 

Hybrid 3.820 4.343 4.749 

Hybrid* 4.010 4.464 4.884 



K = and K = 3 to study this possible effect. As our value of 0:5 is small, the two spectra 
do not differ strongly. 

The colour wave function of two gluons is symmetric and the spin wave function con- 
sidered here is symmetric since 5 = 2. As gluons are bosons, the spatial wave function is 
also symmetric and L is even. The first three states are computed with the six lowest light 
mesons of the first Regge trajectory. Results are given in Table IIIII and in Fig. ^ When 
ft = 3, the trajectory is 

J = 0.36M2 + 0.80. (12) 

When K = 0, the slope is not changed and the intercept becomes 0.57. This results is 
somewhat different from the results of Refs. 0, |2^ because the smaller value of our string 
tension a. But our glueball masses are compatible with those of the model III. A of Ref. p|. 
Moreover, the lightest state has a mass compatible with possible experimental candidates 

IB- 



V. HYBRIDS 

With formula we find k = 3/2 for a gluon-quark pair inside the hybrid. As the quark 
and the antiquark are assumed to be very close, the total colour factor for the Coulomb-like 
force between the gluon and the point-like cc pair is 2 x 3/2 = 3, which is the same as into 
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a two gluon glueball. As in the previous case, we have considered both k, = and k = 3. 

The situation is different for the confinement which is actually a three-body force. Let us 
consider the general case of a triangle formed by three colour sources named A, B, C, and 
a point P inside this triangle, where the three flux tubes, with lengths \AP\, \BP\, \CP\, 
generated by these sources meet jsil]. The triangular inequalities imply \AB\ + \AC\ < 
2\AP\ + \BP\ + \CP\. Let us assume that the energy of these flux tubes are respectively 
o"a|^-P|, cr|-BP|, and cr|CP|, with = ka. The total energy of the flux tubes is then 
(T(A;|y4P| + \BP\ + |CP|). If A; > 2, the above inequality shows that it is energetically 
favourable for the junction point P to be on the source A. The potential energy of the 
system is then cr(|y4P| + |v4C|). In the case of a hybrid gluelump, the gluon is at position 
A with k = 9/4 and the quarks occupy positions B and C, which are merged here. So the 
confining energy is equal to two times the confining energy in a meson, and we take a = 2. 

To compute the mass of the hybrid, we must fix the mass of the cc pair. In our model, 
there is no direct confining interaction between the quark and the antiquark, but a small 
repulsive Coulomb force exists since formula (jHl) implies that n = —1/6 for a quark- antiquark 
pair in a colour octet. The c and c quarks does obviously not really occupy the same position. 
If we assume a typical separation around the J/ip radius, the contribution of the Coulomb 
interaction is around 10-20 MeV. For heavy quark, the constituent mass /i is around the 
current mass m. For a cc, the contribution of the self-energy can be estimated around 
—40 MeV. So, if we neglect the kinetic energy of the quarks assumed fixed at the centre 
of mass, we can take the point-like cc with a mass equal to 2 rric. We estimate that this 
approximation can generate an error on hybrid glueball around 100 MeV. 

The masses of hybrids are given up to J = 4 with the masses of cn and cc mesons (see 
Table Hm and Fig.|2)). The hybrid masses found are significantly lower than those obtained in 
Ref. again because the smaller value of our string tension a. If the Coulomb interaction 
is cancelled, the masses increase by about 150 MeV. 



VI. CONCLUDING REMARKS 



We compute the masses of the lightest glueballs and hybrid gluelumps (gluon with a 
point-like cc pair) within the framework of a simple effective model of QCD derived from 
the rotating string model ^ 3|- The corresponding semirelativistic Hamiltonian is dom- 
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TABLE IV: Values of effective string tension a and colour factor k for different hadrons. 

Hadron cr k 



Meson a 4/3 
Glueball 9a/4 3 
Hybrid 2a 3 



inated by a linear confinement, supplemented by a one gluon exchange interaction and a 
contribution from self-energy for the quarks only. It is assumed that the Casimir scaling hy- 
pothesis is well verified, and that two gluons are attached by an adjoint string in a glueball, 
while the gluon and the colour octet cc pair are attached by two fundamental strings in a 
hybrid gluelump. Values of effective string tension a and colour factor k for the different 
hadrons considered here are gathered in Table IIVI 

In order to minimize the contributions of the spin dependent interactions, only the masses 
of states with maximal spin have been calculated. This strongly constraints the values of our 
three free parameters: the mass of the charm quark, the strong coupling constant and the 
string tension in a meson. With these parameters, we find glueballs and hybrid gluelumps 



with masses relatively small |9|1. The contribution of the Coulomb interaction is about 
100-200 MeV, depending on the particular hadron. The ground state glueballs lie on a 
Regge trajectory and the lightest 2^^ state has a mass compatible with some experimental 
candidates P, Despite the simplicity of the model, these predictions can serve as a guide 
for experimental searches of exotic hadrons. 
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